Chaotic Repellers in Antiferromagnetic Ising Model 



N. S. Ananikian/") S. K. Dallakian, N. Sh. Izmailian, 

and 

K. A. Oganessyan 
Department of Theoretical Physics, Yerevan Physics Institute, Ahkhanian Br. 2, 

375036 Yerevan, Armenia 

February 1, 2008 

Abstract 

For the first time we present the consideration of the antiferromagnetic Ising model 
in case of fully developed chaos and obtain the exact connection between this model 
and chaotic repellers. We describe the chaotic properties of this statistical mechanical 
system via the invariants characterizing a fractal set and show that in chaotic region it 
displays phase transition at positive "temperature" (3c = 0.89. We obtain the density of 
the invariant measure on the chaotic repeller. 
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It is now widely recognized that fractal objects play a fundamental role in many branches 
of physics H, ^. Much progress has been made in finding an appropriate description, in 
terms of "chaotic repellers" or "termodynamic formalism" for these structures (see ref. [^] 
and references therein). The main goal of this Letter is the investigation of the real statistical 
mechanical system in case of fully developed chaos. 

Recently the three-site antiferromagnetic Ising spin model (TSAI) on Husimi tree have 
been investigated. First it was shown that this approach yields good approximation for the 
ferromagnetic phase diagrams , which closely match the exact results obtained on a Kagome 
lattice 0. Second, the quantitative picture of full doubling bifurcations diagram including 
chaos for the magnetization of this system with finite coordination number was obtained 0. 
Note also, that in the anisotropic-next-nearest-neighbor Ising model on a Cayley tree in the 
infinite coordination limit the existence of chaotic phases associated with strange attractors 
have been obtained 

With "termodynamic formalism" we investigate this statistical system in chaotic region 
and describe the chaotic properties of it via the invariants characterizing a fractal set (e.g. a 
strange attractor). In the present work we obtain the strong connection between TSAI system 
in case of full chaos and chaotic repellers. In particular we obtain the density of the invariant 
measure on the chaotic repeller and focused on whether the termodynamical quantities have 
phase transition It is, in general, hard to determine such behavior unambiguously by 
numerical methods if one dose not have further arguments or exact solutions (as in case of 
X —>■ 4x{l ~ x)) to depend on. However, there is a possibility to remedies these deficiencies if 
consider the characteristic Lyapunov exponents as an order parameter, which will differ in two 
phases. Here we manage to calculate this order parameter as a function of the "temperature" 
and clearly show that the phase transition occurs at positive "temperature", whereas in the 
fractal sets, generated by a class of maps close to x ^ 4x(l — x) the phase transition occurs at 
negative "temperatures". 



The pure Husimi tree [jl0[ is characterized the 7-the number of triangles, which go out of 
each site. The Oth-generation is a single central triangle. 

The TSAI model in the magnetic field defined by the Hamiltonian 

A i 

where cXj takes values ±1, the first sum goes over all triangular faces of the Husimi tree and 



the second over all sites. Besides we denote J3 = (3J'^, h = f3h' , (3 = 1/kT, where h-external 
magnetic field, T-temperature of the system and J3 < corresponds to the antiferromagnetic 
case (in all our numerical calculations we will put J3 = —1). 

The advantage of the Husimi tree introduced is that for the models formulated on it, exact 
recursion relation can be derived 0, pJ| : 



— -f( \ -f( \ — ^/^ ^ ~^ -|- z 

where z = e^'^^, /i = e^^ and < < 1. The function f{x) is unimodal: continuous, 
continuously differentiable and has one maximum at x* = fi^^^^^"''^ in [0, 1]. Note that f{x*) = 1 
for any 7, h and T. This function is nonhyperbolic (hyperbolicity for ID maps means that 
1 < /' < 00 in all points) and maps the interval [0, 1] onto [z, 1]. 

The Xn has no direct physical meaning but through it one can express the magnetization of 
the central base site: 

and other thermodynamic parameters, since we can say that the x„ determine the states of the 
system. 

How to link the TSAI system in case of fully developed chaos with strange invariant objects 
(Contor sets) of the phase space, called chaotic repellers or semiattractors? This is the question 
we shall focus on. For this we need a natural partition and that is provided by the cylinders 
(see ref. §]). 

The recursion relation given by eq.(|l|) when 7 = 4 will have the following form: 

X - f(x A fix) - ^^^^^ ^ ^^^^ ^ ^ (2) 

Xn-nXn-l), /W-^23.6 + 2z/Xx3 + l- ^ ^ 

In the case when equation f{x) — x = has only one solution there are no points of sequence 
Xn which can be escaped and therefore no cylinders rise. It means that there is no possibility to 
study the chaotic region of TSAI system in terms of "termodynamic formalism". However in 
some special cases, when equation /(x) — x = have two or three solutions there is a possibility 
to avoid this problem. 

Let us consider the more simple case when equation f{x) — x = have two solutions: xq 
and Xi. So the point xq must satisfy to condition: /'(xq) = 1 . 
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On the other hand we demand that the values of the function eq.(^ coincide at points x = 1 
and X = Xq- From these conditions after some calculations we obtain: 

2 



Xq 



3 + v^- V15 + 6v^ 



(3) 



and corresponding values of fi and z are: 

fx = fiQ = 215.5128, z = zo = 0.01854. (4) 

The second solution of equation f{x) — x = is: 

a;i = 0.3129, = 2.68. (5) 

The plot of function with these values of fio and zq is shown in Fig.l. One can see from Fig.l 
that in this case sequence Xn maps the interval [xq, 1] onto [xq, 1] and the points which come 
to X* = y/xo are escaped. 

As known the analytical criterion for local unstability is: 

One can see from eq.(§) that the mentioned condition of unstability is fulfilled. 

On the other hand following to ref. for 7 = 4 and Zq = 0.01854 we calculate the 
Feingenbaum exponent 6 {6 = 4.669) and also /ioo = 3.18 from which the chaotic behavior 
ensues and sequence essentially never repeats itself. It shows that the value of /xq of eq.@ 
corresponds to chaotic region. 

Thus we obtain the exact connection between the Ising spin model with multisite antifer- 
romagnetic interaction on Husimi tree in external magnetic field (when 7 = 4) in case of "fully 
developed chaos" and chaotic repellers. 

It is interesting to note, that if one lets 7 = 3 in eq.(|I]) rather then 7 = 4, the above 
situation changes dramatically. In this case equation f{x) — x = have only one solution 
in chaotic region. It means, that there is no possibility to study the TSAI system in case of 
"crisis" in terms of chaotic repellers. 

For a crisis map given by eq.(|^) and eq.(§) we want to describe the scaling properties of 
the attracting set, which in this case is the entire interval /: [xq, 1]. For this we need a natural 
partition and that is provided by the cylinders (we follow here refs. [T^). For an index n, I 



is partitioned into 2" intervals or cylinders, these being the segments where the sets of points 
having identical symbolic-dynamics sequences of length n with respect to the maximum point 
X* = ^/xq- "^^^ inverse function h = has two branches and the n- cylinders are all the 
nth-order preimages of /. The length of cylinder is denoted by lei,e2,...e„ , where e G { — 1, 1}. 
The partition function Z{P) is defined as 

Zn{f3)= E e-^'"'^i---, (6) 

€1,62, ■■■en 

where f5 G (— oo, cxd) is a free parameter, the inverse "temperature". In the limit n —>■ oo the 
sum behaves as 

= e-"^(^), (7) 

which defines the free energy, F{P) (the topological pressure JlSl). The entropy 5'(A) is the 
Legendre transform S{X) = ~F{f3) + XjS, where the relation between A and (3 is found from 
A = F'{P), P{X) = S'{\), and these have the following meaning: in the limit n — > oo, e""^^'^-' is 
the number of cylinders with length / = e^'"'^ or, equivalently, with Lyapunov exponent A. The 
dimension of the set of points in / having Lyapunov exponent A is 5'(A)/A. 

With using the eqs.(^, (|^) we can numerically calculate the free energy, which is shown in 
Fig. 2. One can see from Fig.2 that the free energy have a singularity around /? = 1 and shows 
the existence of the phase transition of first order in this region of /3. 

How can the transition be determined accurately? 

Let us consider the characteristic Lyapunov exponent as an order parameter which will 
differ in two phases. Fig. 3 shows this order parameter at positive temperatures calculated for 
different sizes of the system, corresponding to n = 9, 11, 13. The curves converge towards the 
line and result in a first order transition at positive jSc = 0.89, whereas in the fractal sets, 
generated by a class of maps close to x — 4x(l — x) the phase transitions of first order occur 
at negative Pc IJl^ • 



To consider the above obtained results in terms of entropy function S{\), let us first discuss 
the general appearance of the entropy function. First of all it should be positive on some interval 
[Amin, Amax]- The value A = ln2 must belong to that interval, which follows from the fact that 
the sum of the lengths of all cylinders on a given level is 1. Secondly it is often found that the 
values of Xmin and Xmax are given by the logarithms of the slopes at the origin p. In our case 
Xmin = (/'(xo) = 1) and Xmax = ln2.68 (|/'(xi)| = 2.68). Consequently the slope at the fixed 
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point away from xq (xi) is larger then the slope at the origin xq. It also indicates the quite 
different behavior of TSAI system then in the cases of a class of maps close to x ^ Ax{l — x), 
where the phase transition exists ||12|| . 

The precise form of the entropy function is, as mentioned in the introduction, not easy to 
assess with great accuracy. The existence of a first order phase transition implies that there 
should be a straight line segment in S{X). Combining the value of jSc = 0.89 with the point 
corresponding to this jSc'- A = 0.369, ^(A) = 0.414 we find the straight line segment in S'(A) 
when n = 12. Of course, with finite-size data, it is impossible determine the full straight line 
segment in S'(A) and with increasing the n, the straight line segment will increase. 



Finally let us consider the recursion [14 



xe/-i(x') 



In cases with a chaotic attractor, eq.(^) for /3 = 1 (It has been claimed |T3[ that for ID 



maps with chaotic attractors F{1) = 0) is the well known Frobenius- Perron equation ||16|| : 

where Q{x) = lim is the invariant density on chaotic repeller or chaotic attractor. The 

iterative solution of eq.(^ converges very fast {n = 5) and we determine the stationary solution 
of it, which is shown in Fig. 4. 

In summary we have obtained the exact connection between the Ising spin model with 
multisite antiferromagnetic interaction on Husimi tree in external magnetic field (when 7 = 4) 
in case of "fully developed chaos" and chaotic repellers. This remarkable result gives possibility 
to investigate the chaotic properties of this statistical mechanical system via the invariants 
characterizing a fractal set (e.g. a strange attractor). In particular we numerically solved 
Frobenius-Perron recursion equation and obtained the density of the invariant measure on the 
repeller. It is shown that this system in chaotic region displays a first order phase transition at 
positive "temperature" /5c = 0.89, whereas in the fractal sets, generated by a class of maps close 
to X — > 4x(l — x) the phase transitions of first order occur at negative "temperatures". It is, 
in general, hard to determine such behavior unambiguously by numerical methods if one dose 
not have further arguments or exact solutions to depend on. However, there is a possibility 
to remedies these deficiencies if consider the characteristic Lyapunov exponents as an order 
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parameter, which will differ in two phases. In this Letter we manage to calculate this order 
parameter as a function of the "temperature" and clearly show where the phase transition 
of first order will occur, which describe transitions in the distribution of the characteristic 
Lyapunov exponents. This phase transition in terms of free energy and entropy function is also 
analysed. 

The investigation of chaotic statistical physical system has opened new challenges for the- 
ories of stochastic processes, especially in the direction of stochasticity of vacuum in QCD. 
Note that the Z{3) gauge model with a double plaquette representation of the action on the 
generalized Bethe lattice is dual to the spin-one Blume- Emery-Griffiths model |T7| , p!8|] . This 
allows us to suppose that TSAI on Husimi tree when 7 = 4 can be connected via duality with 
multiplaquette representation of the gauge theory. Thus we hope that the present work will 
stimulate the further investigation of stochastic processes. 
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Figure Captions 
Fig.l. The function eq.(|]) for values of /iq and Zq of eq.(^. 

Fig. 2. The free energy F{f3), which have a singularity around (3 = 1 and shows the existence 
of the phase transition of first order in this region of (3. 

Fig. 3. The order parameter X{(3) calculated for different sizes of the system, corresponding 
to n = 9, 11, 13. The curves converge towards the line and result in a first order transition at 
pc = 0.89. 

Fig. 4. Invariant density on chaotic repeller. 
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